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ABSTRACT 


A study to determine the effects of pin elasticity, friction, and 
clearance on the stresses near the hole in a pin-loaded orthotropic 
plate is described. The problem is modeled as a contact elasticity 
problem using complex variable theory, the pin and the plate being two 
elastic bodies interacting through contact. This modeling is in con- 
trast to previous works, by other investigators, which have assumed that 
the pin is rigid or that it exerts a known cosinusoidal radial traction 
on the hole boundary. Neither of these approaches explicitly involves a 
pin. A collocation procedure and iteration were used to obtain numeri- 
cal results for a variety of plate and pin elastic properties and vari- 
ous levels of friction and clearance. Collocation was used to enforce 
the boundary conditions at a finite number of points around the hole 
boundary and iteration was used to find the contact and no-slip regions 
on the boundary. Details of the numerical scheme are discussed. The 
study shows that pin elasticity is not as important as clearance, fric- 
tion, or the elastic properties of the plate in determining contact 
stresses. 


ii 
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Introduction 

Because of the widespread usage of fiber-reinforced composite mate- 
rials in structures, there have been a number of studies to determine 
the stress distribution around a hole in a pin-loaded orthotropic plate 
[1-16], The studies have been two dimensional plane stress analyses 
aimed at understanding the behavior of pinned and bolted connectors made 
of fiber-reinforced composite materials. While the plate has been mod- 
eled various ways, including finite-element [1,5-8,10,12,14-16] and 
elasticity [2-4,9,11,13] approaches, none of the studies have directly 
addressed the pin, or more importantly, pin/hole interaction. For the 
most part the investigations have assumed that either the pin was per- 
fectly rigid [1-9,13-16], or that the pin produced a known cosinusoidal 
radial traction on the hole boundary [11,12,16]. This latter assumption 
was first used by Bickley [17] in his study of isotropic plates. Al- 
though this has been shown to be a good approximation for isotropic 
plates [18], its applicability to orthotropic plates has not been 
checked. Neither approach requires an explicit pin model and each an- 
alysis reduces to a boundary value problem involving a single body, 
namely the plate. The purpose of the work reported on here was to ex- 
plicitly model the pin and its interaction with the hole when determin- 
ing the stresses around a hole in a pin-loaded orthotropic plate. The 
pin model includes pin elasticity, and friction and clearance between 
the pin and the hole. The problem involves two elastic bodies in con- 
tact and is much more difficult than the single body problem. One major 
difficulty is that the regions of contact and no contact between the pin 
and hole are unknown a priori and must be solved for as part of the an- 
alysis. In addition, the presence of friction complicates the problem 
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because the slip and no-slip regions must also be determined as part of 
the analysis. While pin elasticity [10], friction [2,3,5,7,13,14], and 
clearance [2,3,5,15] have been included in previous analyses of this 
problem, no analysis has examined all three simultaneously. The one an- 
alysis which did include pin elasticity used finite elements to repre- 
sent the pin and the plate, the pin transmitting forces to the plate 
through springs connecting the nodes of the pin and plate. In that 
study the effects of pin elasticity were never established. 

The analysis here is cased on formulating the elasticity problem in 
terms of complex variobl° theory. Both the pin and the plate are line- 
arly elastic and the plate is infinite in extent. The pin loads the 
plate through a body force acting on the pin. A solution to the problem 
is obtained by a numerical scheme, namely a collocation procedure and 
iteration. The collocation procedure enforces interface and boundary 
conditions at the pin/hole boundary at a finite number of points. The 
iteration procedure is used to find the contact and no-slip regions. 

This paper begins by formally stating the problem. The key steps in the 
method of solution, which depends on having elasticity solutions for the 
plate and pin to somewhat arbitrary boundary tractions, are then out- 
lined. Following that the major components in the elasticity solutions 
are presented. Next the details of the collocation and iteration pro- 
cedures are described. The number of collocation points, the iteration 
technique, the number of iterations, and solution accuracy are di- 
scussed. Finally, numerical results are presented. These results illus- 
trate the effect of pin elasticity on the stresses around the hole and 

also indicate the importance of pin/ hole friction and clearance on the 
stress calculations. The effects of the plate's elastic properties on 
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the stresses are illustrated with the computations. The paper closes 
with a brief discussion. 

Statement of the Problem 

The plate's principal material axes are aligned with an x-y coordi- 
nate system. The origin of the x-y system is the center of the hole, as 
is the origin of a cylindrical r-9 coordinate system used in the analy- 
sis. The hole is of unit radius and the plate is of unit thickness. The 
pin loading the hole is isotropic and the net force the pin exerts on 
the hole is in the x-direction. The pin radius is 1 - a and, due to 
loading the hole, the center of the pin moves a distance 6. Figure la 
depicts the geometry and coordinate systems used in the analysis. Fig- 
ure lb illustrates the parameters X and 6. The case of X equal to zero 
represents a perfectly fitting pin. As will be seen, this is a special 
case. A +6 is in the direction of the +x axis. Here <5 is considered 
the independent loading parameter. 

As the pin moves to the right, it contacts a portion of the hole. 
The extent of the contact region is unknown and the half contact arc is 
denoted as 3. Due to friction within the contact arc there is a region 
of slip and a region of no slip. The half no-slip arc is denoted as a. 
Outside the region of contact there is a no-contact region. Figure 2 
illustrates the various regions on the pin/hole boundary. Coulomb fric- 
tion is assumed to act in the contact region. The coefficient of fric- 
tion is denoted as u. It is assumed a priori that there is one contact 
region and within that region there is one no-slip region. It is con- 
ceivable that there might be multiple regions of contact and no-contact, 
or slip and no-slip. However, when the stresses are actually computed 
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tr.j results can be examined for this possibility. If there are multiple 
regions of contact, tensile radial stresses will most likely develop in 
those regions where contact is erroneously assumed. For the cases con- 
sidered there has been no evidence of this sort of behavior in the re- 
sults. Also, examining the friction-induced shear stresses and their 
relation to the radial stresses will indicate whether multiple no-slip 
regions might occur. There has been no evidence of this behavior 
either. It is assumed that the coefficient of friction is independent 
of circumferential location. This assumption may not be as valid for 
composite materials as it would be, say, for metal. 

Formally the boundary and interface conditions between the plate 
and the pin can be stated as follows: 

In region I, the no-slip region, -a < 9 < a: 

u p = u r + 6 cos 9 - \ (1) 

r plate pin 

u = u - 6 sin e . (2) 

w plate 9 pin 


In region II, the slip regions, a < 9 < p and - p < 9 < - a: 


u = u r +6 cos 9 - \ (3) 

plate pin 


lv 9 l = 


In region III, the no-contact region, p < 9 < 2n - p: 


°r = 0 = Ve* ( 5 )» 

r re (6) 
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In the above the radial and tangential elastic displacements of the 
pin and plate are given, respectively, by 


pm 


pin 


plate plate 


Equations 1 and 3 express continuity of radial displacements in the 
contact region and eq. 2 expresses continuity of tangential displace- 
ments in the no-slip region. Equation 4 is the Coulomb friction law. 
Equations 5 and 6 express the traction-free conditions in the no-contact 
region. Because of the nature of friction, assumptions have to be made 
regarding the sign of x r0 . The problem can be solved assuming 6 is 
increasing. This corresponds to a positive t^. Alternatively, the 
problem can be solved assuming 5 is decreasing. This corresponds to a 
negative c rlJ . Only the former condition will be discussed here. 


Method of Soluti on 

Due to equilibrium considerations the boundary tractions acting on 
the pin are identical to the boundary tractions acting on the hole. As 
formulated here the problem reduces to finding these boundary tractions. 
Once the tractions are known the stresses in the plate can be deter- 
mined. Here the unknown boundary tractions are represented in the form 
of a complex Fourier series with unknown coefficients A^. The series is 
given by 


(N-n) pin = (H- 


plate 


+oo 

V 

L 

-00 


iko 

e 


all 0 


(7) 


where, following the convention of [19], 
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N = the normal, or radial, traction a r on the boundary, and 

T 3 tangential, or shear, traction r . on the boundary. 

ry 

Because the tangential tractions are odd functions of 6 and the normal 
tractions are even functions of 9, the A k are real. 

i kG 

Since each is a legitimate t action, each term e in the series 

produces a unique set of stresses and displacements in the plate and in 

the pin. In the next sections the stresses and displacements in the pin 

i k0 

and plate due to traction e will be found. If the boundary traction 
i kG 

is A k e , then the stresses and displacements in the plate and pin 

will be weighted by A k . If, for example, the radial displacement due to 
i k G 

N-i T = e is denoted as (u p ) k , then the radial displacement due to 

i kG 

traction N-i T = A k e is given by A k (u p ) k . Considering all terms in 

eq. 7, the radial displacement u r would be the weighted sum of the ef- 

i ko 

fects of the individual traction? c . That is. 


u 


r 



V u r>k 


( 8 ) 


If all A k were known, then u p could be calculated. The same can be said 
of the tangential displacement, namely 


% % J Wk > (,) 

K = -<2> 

iko 

where (u Q ) k is the tangential displacement due to traction N-i T = e 

No distinction has yet been made between the plate and the pin. The re- 

i k 0 

sponses of the plate due to N-iT = e are different than the responses 

i k G 

of the pin due to N-iT = e 


Thus the notation 




- I \(u r ) k 

k=-» * r pin 


7 
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Wk 


pin 


( 10 ), 

( 11 ) 


plate 


= I \( ,j r ) k ; 

k=-« K r K plate 


'9 


plate k=-® 


l Wk 


plate 


( 12 ), 

(13) 


distinguishes between pin and plate responses. It is important to point 

i l/Q 

out that since they are a result of traction e , the quantities 


< U A 


pin 


<Vk 


pin 


<«A 


plate 


and (uj 


e'k 


plate 


are assumed to be known. The problem will reduce to finding A k . 

To continue with the formulation of the problem, the stresses in 
the pin and plate are given by 




plate 


\l ** ( v- 


plate 


(14), 

(15) 



k 


2 A k 


pin 


plate 


k 


l \<Vk 


plate 


(16), 

(17) 



k 


i y 


T re^k 


pin 


're 


plate 


• l 


V T re*k 


plate 


(18), 

(19) 


The quantities (a). (t ). are the stresses in the pin 

pin ik0 plate 

and plate due to traction N-iT = e . They are also considered known. 

To determine the unknown coefficients A^, a collocation technique 
and iteration are used. Instead of using infinite sums in eqs. 10-19, 
finite sums are used. Instead of satisfying the boundary conditions 
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given by eqs. 1-6 at each point oo the circular boundary, (i.e., at an 
infinite number of points) the boundary conditions are satisfied at a 
finite number of points. The number of terms in the truncated series is 
relaced to the number of collocation points on the boundary. The itera- 
tive portion of the solution involves finding the contact and no-slip 
arcs, a and p. They, like the A k , are unknown. However, if specific 
values are assumed for a and p, a solution to the problem with those 
values of a and p can be found. The specific values may not yield a 
solution which satisfies all the boundary and interface conditions, but 
iteration can be used to find values which do satisfy all the conditions 
of the problem. The iteration procedure will be further discussed. 
First, the collocation procedure is explained. 

By using finite sums in eqs. 10-19 and substituting those sums into 
the boundary conditions for the various regions, eqs. 1-6, the method 
for determining the A k 's becomes evident. Substituting the finite sums 
into eqs. 1-6 results in: 

For the no-slip region; 


N 

I \( u r )k 

k=-N * r K plate 


l A k (u ) 

k=-N K 9 K plate 


N 

l \(u p ) k + 6 cos e - x 
k*-N pis 

N 


■ l 

k=-N 


Wk 


pin 


- 6 sin 0 . 


(20) 

( 21 ) 


For the slip regions; 


N 


l 

k--N 




plate 


N 

3 )’ M U A + 6 cos 0 - x 
k=-N K pin 


k J N A k^ x r0*k " *■ 0 * 


( 22 ) 


(23) 
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For the no contact region; 

N 


I A k (oJ k = 0 , 
M 

(24) 

N 

J M A kWe>k = 0 * 

(2b) 


There are 2N+1 terms in each sum. If the boundary conditions are satis- 
fied at 2N+1 discrete points around the circular boundary, a set of 2N+1 
linear algebraic equations can be formed from which to solve for the A k . 
Since the proper symmetry conditions have been incorporated into the 
solution, the 2N+1 discrete points need only be distributed over one- 
half the hole circumference, 0 < 0 < u. The assumed values of a and p 
determine over what range of 9 eqs. 20 and 21 are enforced, over what 
other range eqs. 22 and 23 are enforced, and over what range eqs. 24 and 
25 are enforced. Since the tractions on the plate and pin are identi- 
cal, neither the plate nor the pin is called out specifically in eqs. 24 
and 25. Though it is not specifically noted, attention must be given to 
signs in eq. 23. In addition, strictly speaking, by eq. 7 


<*r>l 


pin 


(oj 


r k 


plate 


= cos(ke) 


and 




pin 




plate 


= - sin(ko) 


(26) 

at the 
boundary 

(27) 


Equations 23-25 thus can be written as 
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N 


l A k (slu(ke) + p cos(k©)) = 0 . 

(23)a 

-N 


N 


l A. cos(ke) = 0 
k=-N 

(24)a 

N 


l A. sin(ke) = 0 . 

i* hi K 

(25)a 


Expanding eqs. 20-22 and 23a-25a, evaluating them at 2N+i circumferen- 
tial (e) locations, rearranging slightly, and putting them into matrix 
notation ieads to equations of the form 


[C] {A} = {B}. (28)a 


Expanding, these equations can be written as 



cos(e)- 
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The terms in the top two lines in the coefficient matrix [C] and in 
the vector (3} represent terms evaluated at various 0 locations in the 
no-slip zone, the terms in the second two lines represent terms evalua- 
ted at various e locations in the slip zone, and the last two lines rep- 
resent terms evaluated at variojs 0 locations in the no-contact zone. 

It is important to point out that with a perfectly fitting pin \ = 0 and 
each element on the rignt-hand-side of eq. 28 changes, for example, by a 
factor of two when 6 is doubled. Each A^ simply doubles and the spatial 
variatio' of N-iT, with 0, remains unchanged. Thus a and p are unaf- 
fected by c. qith x t 0, doubling h does not result in a doubling of 
the right-hand-side of eq. 28. Each A k changes by a different amount 
and thus the spatial variation of N-iT varies. This translates into 6 
directly affecting a and p. The perfect fit case is therefore quite 
special and is often referred to as the linear case. 

To implement the solution procedure the responses of the plate and 

j |^0 

pin to boundary traction e , namely (u ). , (u ). must 

r K pin K plate 

be found. This follows. 

i ke 

Res pons e of Pin to e 

In reality, the pin transmitts a load to the hole boundary because 
it is loaded by some external mechanism. Commonly the pin is in double 
shear, being loaded by two other plates parallel to and on either side 
of the plate under discussion. This is shown in fig. 3a. To be sure, 
there are bending stresses in the pin and the exact determination of 
stresses in the pin requires a three-dimensional analysis. However, for 
the segment of pin within the thickness of the plate the situation shown 
in fig. 3b is quite accurate. This segment of the pin, shown as a disk. 
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is loaded by shear stresses from portions of the pin on eitner side of 
the plate, and by the contact stresses. Compared to the length of the 
pin, region of the pin within the thickness of the plate is small. 
However, the fact that the pin actual ij :.u; ’“noth does influence the 
deformation of the pin within the thickness region. Here the pin is as- 
sumed to be in a state of plane deformation and the shear stresses in 
fig. 3b are assumed to be spatially uniform over the pin cross-sectional 
area. Furthermore, for purposes of analysis, the effects of the uniform 
shear traction are assumed to be equivalent to a uniformly distributed 
body force acting in the x direction within the pin. An alternative pin 
model would be to assume the disk representing the pin is loaded at its 
center by a concentrated force. This model was not considered because 
it was felt to be too harsh an idealization of the forces actually 
transmitted to the portion of the pin within the thickness of the plate. 
A concentrated force model was adapted by Rao [20] in a study of iso- 
tropic plates. 

With the model adopted here, the equilibrium equations for the pin 


are 

d°r 1 0t r 9 

dr r 09 


+ b„ = 0 , 4 


OO, 


00 


Ox 


ro 


2x 


re 


dr 


b n = °» 
w (29)a,(29)b 


where b and b are the components of the body force. Mi 1 ne-Thompson 
[21] discusses the complex variable approach for the case of a body 
force and so only the primary features of the analysis are presented 


here. 


To effectively use complex variables, it is convenient to express 
the body force components in terms of a scalar potential, V, i.e. 


b 

r 


= _ 0V(r,9) 
dr 


and b„ 
9 


1 dV( r ,Q) 

r on 


(30) 
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Defining a stress function F(r,e) as 


o 


r 


I &£. + I_ + 

r dr r 2 d6 2 


V; 



(31)a,b 


* l_ a 2 _F_ _ 1 d 2 _F 
r9 r 2 ae 2 r arae 


(31)c 


automatically satisfies the equilibrium equations. Using Hooke's Law 
for an isotropic material, with Poisson's ratio v, and substituting the 
strains in terms of the stresses into the only nontrivial compatibility 
equation leads to the following equation for F and V: 


where 


Defining 


V 4 F + kV 2 V = 0, 

K - (1 - 2v)/(l - v). 
V * V 2 Q, 


(32) 

(33) 

(34) 


where Q is another scalar function, eq. 32 becomes 

(35) 

(36) 


Defining 


V 4 (F + kQ) * 0 . 


ie 

i = x + iy = re , 


the solution for F can be expressed in terms of two complex functions, 
$ and x. and Q as 
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F = | [z*(z) + zl('z) + x (z) + x(z)] - tcQ(z,z) . (37) 

The first and second fundamental stress combinations are written as 

o + o r = 2($(z) + <l>(z)) + (8 - 4 k) , (38) 

dzaz 

a „ - <v + i?x fi « 2(z$‘(z) + y(z) - 2 k — e 2 ’ 8 , (39) 

J ' 8 5Z2 

and the elastic displacements are given as 

( u r + l ’ J e ) = (O - 4v))*(z) - zt'(z) - l(z) + 2k -- ^ 2 ^ )e~ 16 . (40) 

dz 

In the above the prime, denotes differentiation with respect to the 
argument, G is the material's shear modulus, and 

$ z) * <p‘(z) , f(z) = 4 ,'(z) , and ,p(z) = x ‘(z)* (41)a,b,c 

The overbar, , denotes complex conjugate. 

At the boundary 

"r * ’Ve * " ' ” \ I V'*'- («) 

Therefore the stress cou. ination most useful for matching traction 
boundary conditior is obtained hy subtracting eq. 39 from eq. 38. Tne 
result is 

°r ‘ lT re = + " (**'!*) + y(z))e 2ie 

♦ i(2 - „) sMlJl t 3.11UJ1 e f2e_ 

hzoz a 2 z 


( 43 ) 



16 


orig^al r :-. r r : *3 

OF POOR QUriwlTY 

Since the pin domain is simply connected, $(z) and Y(z) can be repre- 
sent en by 

«(z) = l a z k - l a k r k e 1k0 (44)a 

k=o k=o 

*(z) = l b .z k « I b.r k e ik0 , (44)b 

k=0 K k=0 K 

a k and b k being constants, possibly complex. 

To complete the pin analysis, Q needs to be determined. Milne- 
Thompson [21] presents a derivation of Q for a more general body force 
and those results can be reduced to the situation here. The result is 

b 

U = - (z 2 z + (z) 2 z)» (45) 


where b„ is the body force in the x direction, b being given by 

X A 


b x = ^ b r + b e • ( 46 ) 


Substituting eqs. 44 and 45 into eqs. 38-40, the stresses and displace- 
ments can be determined in terms of a k , b k and b x . Specifically the 
stress combination of eq. 43 evaluated on the boundary is 


k*« 


“ a ikO _ X q k i ko . T - „k -iko 

). A. e * l a. R e + l a.R e 


k=0 


k=0 


- I ka k R k e ik3 - l b k 2 R k ' 2 e ik0 
k=l K k=2 K * 


b R(2 - k) .. b Rk . 

-i- (e 10 +e - l0 J --*-e'°. 


4 


(47) 
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Here R is the radius of the pin. Since the hole in the plate is of unit 
radius and the clearance between the pin and hole is x, 

R = 1 - X . (48) 


pfcCsfc # 
QUAUfV 


Using eq. 47 and setting all A^'s to zero except one, and setting 

that equal to one, constants a fc and b^ can be found by matching terms 

on the left and right sides of eq. 47. The solution of the pin to trac- 
ike 

tion N -iT = e is now available. Table 1 shows the result of match- 
ing terms. The values of a fc and from Table 1 can be substituted into 
eq. 44 and those results into eq. 40 to find explicit expressions for 

(u ). and (u ). . These expressions, in turn, are used in eq. 28. 

r k 0 k 

pm pm 

It should be mentioned that b x is not an independent constant. The 
body force and the boundary tractions must be in equilibrium. Except 

i k 0 

for k = +1, boundary tractions of the form e dc not produce a net 
force on the boundary. Therefore, b x is identically zero except for the 
case when Aj is nonzero. For the boundary traction A^e 10 , 



(49) 


It should also be mentioned that in determining the results of Table 1, 
the imaginary part of a Q was arbitrarily assumed to be zero. 

i k 0 

Response of Plate to e 

The plate is assumed to be in a state of plane stress with trac- 

N - iT = l A.e lke 
k=-® K 


tions of the form 
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Table 1 


Coefficients of the stress functions $ and 4, for pin. 


nonzero nonzero 

A k a‘s 


nonzero 

b's 


A 


D 


1 


1 

*1 = 1 

A_2 = 1 


A 


2 


1 



a o = 1 


a l ^ R 


a l = 


a 2 " 


1 

k - 2 

2 R 

1 


R 2 

al 1 a's = 0 


a 3 = 


1 


R 3 

all a's = 0 


all b's = 0 
all b's = 0 
all b's = 0 
b 0 = -l 





A +l( = 1 all a's = 0 


k-2 


b k-2 ‘ 


-(A _-_Jl 

R k-2 

-1 

R k-2 
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acting on the hole boundary. These surface tractions are in equilibrium 
with vanishingly small stresses at infinity. The plate's principal ma- 
terial axes are aligned with the x and y axes and so the constitutive 
behavior can be written as 



Here E x and are Young's moduli in the x and y directions, respective- 
ly* v X y is the Poisson's ratio relating contraction in the y direction 
to the extension in the x direction caused by a tensile stress in the x 
direction, and G is the shear modulus in the x-y plane. The equilib- 
rium equations in the principal material system are 


do 5x 

x + = 


dx 


ay 


3x do 

__*y + _j t 
3x dy 


(51) 


Defining a stress function as 


o 


X 



d£F . 
ax 2 


a 2 F 

T xy “ * axay 


(52) 


the equilibrium equations are satisfied. Using the constitutive 
equations and the stress function, tne only nontrivial compatibility 
equation can be used to determine F. This equation is 


ay 4 b xy x<y dx 2 by 2 E y bx 4 


( 53 ) 
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Lekhnitskii [22] and Milne- Thompson [21] discuss the general solution to 
this equation and associated boundary conditions, while dejong [4,9,13], 
Oplinger and Gandhi [2,3], and Mangalgiri [7] discuss the solution in 
the context of a plate with a hole. Therefore, details of the solution 
are omitted here. 

The solution to eq. 53 has the form 

E = 4“ [ ) + ^2 ( Z 2 ^ + ^1 ( Z 1 ^ + ^2 ^ z 2 ^ * ( 54 ) 

where 

z 1 = x + ujy and z 2 = x + p ? y (55) 

and ^ and p^ are two roots of the characteristic equation associated 
with eq. 53, namely 


E x * 

^ + — " 2v xy^ 2 + E~ = 0 • (56) 

xy J y 

(There are four roots to this equation but they occur in complex conju- 
gate pairs, i.e. p^ = pj and p^ = p 2 ) . The functions and are to- 
be-determined functions of the two transformed variables Zj and z 2 . 

With the definitions of eq. 55 


z k = V + 6 k 2 ’ k = 1 » 2 » 

where 

1 - ip k 1 + ip k 

Y k = ~ dnd \ = — 2 * 


(57) 

(58) 


The first and second fundamental stress combinations are 
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crjok-v-*. . ' ; :*j 

OF FOCrv 

°r + °e = Y l 6 l W l( 2 l^ + y 2 6 2 W 2^ z 2^ + Y l^l W l( Z l) + 

y 2 6 2 W^(i 2 ) . f 59) 

and 

°0 ' °r + 2it r0 = ( Y l W i'( 2 l) * y 2 W 2^2^ + H 

+ l\ W"Ci 2 ))e i20 . (60) 

The displacements are 

( u r + iu e) = itPi W :( z i) + P2 W 2^ Z 2^ + VJCV + J 6 ' 19 * (61) 


with 


Pi = 


Po = 


<?i = 


v xy **1 i ,1 2i 

- - + — + — If" " e 4 l) 


Pi 


4 


_m + 'jL + i-,l *X 21 

E.. E„ u„ ( E. E„ *2> 


-x »*2 ~y 

.2 


+ il + L_ rl- . 521 

E x E x „ (e v E x ^ 


^9 


E E - ''E E ^2 

x c x p 2 y y 


J 


(62) 


The stress combination a f - it r(5 can be formed, as was done with 
the pin, to match the traction N - IT at the boundary. Howev<*'\ 
definitions of eq. 55 transform the problem from the z plane to the 
and z 2 planes. These transformations distort the unit circle (hole) of 
the z plane into ellipses in the z^ and z 2 planes. The boundary condi- 
tions are not defined in these planes. A second set of transformations 
which map the z^ and z 2 elliptical boundaries into unit circles is 
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necessary. These transformations map the Zj plane onto the Cj plane and 
the z 2 plane onto the c 2 plane. These transformations are 


r _ Z 1 + / 2 1 - Vi 

2Yl 


Z 2 + ^ Z 2 ~ 4 y 2 6 2 
C 2 ‘ 2y 2 


(63) 


i 9 i 6 

As can be seen, when z = e » Cj = C 2 = e * thus preserving the shape 
of the unit boundary. The inverse transformations are 


z i * m i( c i> ” n c i + ^ 


Z 2 = m 2 ^ 2 ^ = y 2 r ’2 + 


(64) 


With the change of variables, new nomenclature can be introduced, 
namely , 

W^(z k ) = W‘(z k (c k )) 5 20 k (c k ), k = 1,2 . (6b) 

With this 


W 1 *) ■ • k = 1 - 2 


( 66 ) 


The functions 4>| C (Z| < ) are of the form 


°° a. 


* a o’ n 'l ♦ X -J 


C 


k 2 (C 2 ) 3 b 0 ^C 2 + l 


!=2 C 


k ’ 


( 67 ) 
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; *here the a k 's and b^'s are constants. The logrithmic terms reflect the 
fact that there is a net force on the boundary of the singly connected 
plate. 

For the plate the analog to eq. 43, evaluated on the hole boundary, 
is 

(°r" iT re)| e ie s + " ^ Y lW + Y 2 4'2^2^ e ^ I g i 


Using eqs. 42 and 67, eq. 68 becomes 


( 68 ' 


n.ike - io r,- i(k+l)o^ 

l \e - «i(V> -J, k V ) + 

k=-» k=l 


- X - W - ,1 H*- i(k - 1)e ) 


T^,." - l * . 

60 k=l K 


(69) 


As in the procedure for the pin, setting all A k to zero except one gives 

equations for determining a k and b^ and hence the response of the plate 

i kG 

to traction N=iT * e on the hole. However, with the plate the pro- 
cedure is slightly different. No matter which A k is being considered, 
eq. 69 always yields an equation for a Q and b Q , in addition to other 
equations involving other a's and b's. This equation is 


6 l a o + Vo 


5 - • Y l a c 


Y 2 b o 


1 , when A. = 1 

(70) 

0, all other cases. 


More information Is needed to determine a Q and b Q . The information can 
be obtained by examining the displacements. 
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Because of the logrithmic term, there is the issue of single val- 
uedness of the displacements. Using eq. 65 in eq. 61 and evaluating eq. 
61 at 9 = 0 and 0 = 2n, the following equation results: 

K + %>l , = (p l a o * p 2 b o ' Vo • Vo' 2 '* * (“r * 1u e'l „ (71) 

|9=Zit 1 0=0 

For single valuedness, the expression in parenthesis must be zero, i.e. 


Vo + Vo - Vo ‘ Vo ’ 0 


Equations 70 and 72 provide the necessary equations to determine a Q and 
b 0 . Except for k = 1, a Q and b Q are zero. For the case k = 1, eqs. 70 
and 72 were solved for numerically. Table 2 lists the coefficients a k 
and b^ for various values of k. With Table 2 the solution to the plate. 


namel j the responses 
eq. 28. 



aod (u 0 ) k 


plate 


are available for use in 


Oetails of the Numerical Scheme 

The collocation procedure included the points 9=0 and e=u. Because 

of symmetry, the elasticity solutions for the pin and the plate due to 
i k 0 

traction N-iT * e automatically yielded 

(uj k = (u ). = 0 0 0=0, k=U,+l ,...,+N . (73) 

pin plate 

and 


• T ro^k = ^ T r0^k s 0 0 Qsn > k=0,+l,...,+N . 
pin plate 


( 74 ) 
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Table 2 

Coefficients of the stress functions and ^ for the plate 


nonzero 

A., 


nonzero 

a's 


nonzero 

b's 


A o = 1 


A -1 • 1 


A+i " 1 


1 6 2 Y 1 " 6 l y 2 


*2 = 2(6 2Yl - 6^ b 2 = 

soln to eq. 70 and 72 


A _2 * 1- a_ = 


3 3(6 2 y i - 


A + p =1 a = r 

1 6 3 y 1 ‘ 6 1 y 2 


6 1 y 2 “ 

5 2 y 1 

t 

'1 

) 

- 6 2 Y l) 



3l6 lY2 

~ f> 2 Y l") 

Y 1 


6 1 y 2 * 

6 2 y 1 


A -k " 1 


k+1 ‘ (la] 


6 lY 7T k+ i nr+ijtve - 5 2 y i 


, 

5 2 y 1 ‘ 5 1 y 2> 


:-l {k-l)(6^ Y 2 ” ^2 y 1' 
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These two conditions led to two columns and rows of eq. 28 being zero. 
The collocation procedure included points at 6*0 and 6*n because there 
was one other condition at each point that had to be enforced, namely 


and 


( u r)p1h * ( u r)plate " ® at e *° 


a r * 0 0 0 * n . 


(76) 


Thus the two rows and columns representing eqs. 73 and 74 were elimina- 
ted from the set of linear equations, eq. 28. Failure to do this re- 
sulted in a singular matrix. The number of equations and unknown A^'s 
solved for was reduced to 2N-1. 

Values for a, p, and 6 were chosen. The system of 2N-1 equations 
was solved for the A^'s. In all likelihood the values of a and p chosen 
did not satisfy all the interface and boundary conditions for that spe- 
cific value of 6. An iterative process then began to find values of 
a and p that did. The key to the iteration process was the satisfaction 
of conditions which were a combination of a restatement of the condi- 
tions stipulated by eqs. 4-6 and statements regarding the physical 
realities of the problem. These conditions were: 


„ < 0, all e 

(77)a 

. » T r0 - p < e < * 

( 77 )b 

, * T re * 0, 0 * p 

(77)c 

re ! < ^l^pl . 0 < 6 < a 

(77 )d 

r{) l = nk r |. 8*«t. 

(77)e 


at the hole edge 
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Equation 77a states that the radial tractions had to be either compres- 
sive or zero, everywhere. Equation 77b states the tractions had to be 
zero in the no-contact zone. Equation 77c states that at the end of the 
contact zone the tractions had to be zero. Equation 77d states that the 
friction-induced shear traction had to be less than the Coulomb limit in 
the no-slip zone. Equation 77e states that at the end of the no-slip 
zone, the friction induced shear stress was exactly the Coulomb limit. 
Equations 77c and 77e are written explicitly : ecause behavior of a r at 
e * p and t at 9 * « were the most powerful indicators of convergence. 

To obtain a solution, the value of « was set to a value correspond- 
ing to the first collocation point away from 0 = 0. Attention was fo- 
cused on p by selecting a realistic value. Depending on whether the 
value of p chosen was larger or smaller than the correct value, the o p 
vs. 0 relation had very distinctive characteristics near 0 = p. Figure 
4 illustrates the nature of the solution for three values of p. The 
correct value of p is denoted p*. If the value of p chosen was larger 
than p*, say p = Pj, o p was tensile near 9 * p^. If the value of p 
chosen was smaller than p*, say p * p£, then o p was not zero at 0 = p ? . 
By recognizing what the character of the o f vs. 0 meant relative to the 
assumed value of p, within 3 or 4 iterations a value of p could be ob- 
tained which satisfied eqs. 77a-c. Iteration for a then began. This 
iteration procedure was similar in that the x p0 vs. 0 relation showed 
distinct characteristics near 0 * o. These characteristics depended on 
whether the value of a chosen was larger or smaller than the co. ect a. 
The characteristics near 6 = a are shown in fig. 5. The correct value 
of a is denoted as a*. If the value of o chosen was larger than 
a*, say a^, then at 0 * a^, |x | exceeded p|o p j. If the value of a 
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chosen was smaller than «*, say a * « 2 » then at 0 * a^, |x r0 j was less 
than |i|o r | . Within 3 or 4 iterations a solution for a was obtained 
which satisfied eqs. 77d and c. The value of 8 was then rechecked and 
adjusted if necessary. Generally there was little coupling between the 
two variables. 

For smaller values of p (e.g. p = 0.2) the no-slip region was found 
to be small. Often the no-slip region was smaller than the spacing be- 
tween the collocation point at e * 0 and the point next to it. In these 
cases the no-slip condition was satisfied only at 9 * 0. Since by the 
symmetric nature of the solution there was no relative tangential motion 
at e = 0 anyway, the no-slip region could only be determined to be less 
than the collocation point spacing. 

Overall, the convergence of the solution was found to be dependent 
on the location and number of collocation points used to satisfy the 
boundary conditions. It was important to have the collocation points 
close to being equally spaced around the hole. Forty collocation points 
around the half-circle were used to obtain the numerical results pre- 
sented in the next section. Thus the results to be presented used 78 
unknown Aj.'s in the series representing N-iT. The 78 terms were divided 
evenly between negative k and positive k. Twenty to 30 collocation 
points gave reasonable results when compared to similar cases from pre- 
vious investigations. Forty points represented the upper limit of the 
interactive computer used in the analysis. Fewer than 20 points result- 
ed in pour agreement with established work. 
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Numerical Results 

Numerical results were obtained in order to assess the effects of 
the various parameters on the stresses around the hole. The numerical 
study considered plates with varying degrees of orthotropy. The results 
for four different plates are presented here. In the context of fiber- 
reinforced composite plates, these four plates represent four different 
laminates. The material properties of each plate Te indicated in Table 
3, as are the laminates they represent. The plate properties considered 
coincide with the basic properties used by Crews et al. [10] and repre- 
sent a graphite-epoxy material. 

Plate A represents an extreme in degree of orthotropy. The plate 
is much stiffer in the load direction than it is perpendicular to the 
load direction. This represents a graphite-epoxy plate made with all 
the fibers in the direction of the load. Plate 8 represents the other 
extreme, being much softer in the load direction than perpendicular to 
the load direction. Plate B represents a graphite-epoxy plate with all 
the fibers perpendicular to the load direction. Plate C represents a 
plate with properties representative of laminates which have the same 
inplane stiffness in all directions. Such laminates are referred to as 
quasi -isotropic laminates. Plate D represents a moderately orthotropic 
laminate, being three times stiffer in the load direction than it is 
perpendicular to the load. These four plates will be used to illustrate 
the effect of pin flexibility and the plate's elastic properties on the 
stresses around the hole. Plate 0 then will be used to illustrate the 
effect of friction, clearance, and pin displacement on stresses. 
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Table 3 

Material Properties of Plates 


Plate 

E x 


G xy 

V 

representati ve 


Msi 

Msi 

Msi 


laminate* 

A 

21.3 

1.58 

0.930 

0.38 

0° 

3 

1.58 

21.3 

0.930 

0.028 

90° 

r 

w 

8.40 

8.40 

3.20 

0.310 

(07i457‘-0°) s 

D 

12.4 

3.73 

3.21 

0.667 

(0° 2 /t45°) s 


*fiber angles relative to +x direction 

Figures 6-9 show the effects of pin flexibility and plate material 
properties on the stresses around the hole. The three pin flexibilities 
used were: a rigid pin, a steel pin (E * 30 x 10^ psi), and an aluminum 

pin (E = 10 x io 6 psi). Poisson's ratio of the pin had very little ef- 
fect on the numerical results and a value of 0.3 was used. In each of 
figs. 6-9 the coefficient of friction between the plate and the pin was 
assumed to be 0.2. This represents a reasonable value for metal on 
graphite-epoxy. The non-dimensional clearance between the pin and 
plate, a/R, was 0.01 and the nondimensional pin displacement, 6/R, was 
0.03b. Realistic values of pin/hole clearance vary from application to 
application but A/R = 0.01 is representative. Each figure has the same 
scale, for easy comparison, and each figure illustrates the three stres- 
ses at the hole edge, i.e. the radial stress the circumferential 
stress o q , and the friction-induced shear stress T r0 . 


The stresses have 
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been nondi mens ionali zed by the bearing stress and the behavior of the 
stresses with circumferential location is illustrated. The bearing 
stress is defined as the pin load, P, divided by the product of plate 
thickness and hole diameter. Also shown in figs. 6-9 is the cosinusoid- 
al radial stress distribution often assumed. Thus the figures illustrate 
the two assumed pin/hole interaction models used in previous investiga- 
tions, i.e. rigid and cosinusoidal, as well as the current more realis- 
tic one. As a check on the numerical solution, the stresses 
and t pQ were integrated around the hole edge to determine the total 
load acting on the hole. The integration was always within 0.01% of P. 

When examining figs. 6-9, three conclusions are obvious. First, it 
is clear pin flexibility is not a big factor in determining the stresses 
at the hole edge. Until now, no investigation has shown this explicit- 
ly. Second, the degree of orthotropy strongly influences the peak 
stresses and the distribution of the stresses around the hole. The 
highly orthotropic plate A in fig. 6 has a stress concentration factor 
of 2 for the circumferential stress compared to the 1.2 of the quasi- 
isotropic plate C in fig. 8. Third, the often-assumed cosinusoidal dis- 
tribution is not generally accurate. It is a serious misrepresentation 
for plate B and, due to pin/ hole clearance effects, does not properly 
represent the contact region in any situation. It will be seen that the 

character of c near 0 * U is determined by friction while the character 
r 

of o p near 90° is determined by pin/hole clearance. Thus for particular 
friction and clearance levels, the cosinusoidal distribution could be a 
good representation. However, in all cases, the rigid pin assumption is 
better. In all the cases shown in figs. 6-8, the ci rcumferential 
stress a 0 is negative at 0 * 0. This is somewhat counter to intuition 
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but the result, as will be seen shortly. Is due to frictional effects. 
This was found in Wilkinson's [5] analysis and it has been measured 
[23]. The decrease in the magnitude of o r as 0 + 0, In plate A for 
example. Is also due to friction. 

Table 4 indicates the contact and no-slip angle determined by the 
iterative procedure for figs. 6-9 and for figures to be discussed. The 
location of the maximum hoop stress is also Indicated in the table. The 
location of the maximum hoop stress, while generally occurring near the 
end of the contact region, depends on the plate's elastic properties. A 
plate representing a laminate with fibers at +45° and -45° relative to 
the load direction experiences a maximum hoop stress at roughly e s 45°. 

Figure 10 illustrates the effect of friction on the stresses at the 
hole edge of plate D. Three values of friction were examined for the 
case of a steel pin with clearance \/R = 0.01 and displacement 6/R a 
0.035. The values of friction examined were: p » 0, 0.2, and 0.4. The 

1st and 3rd values of p represent extremes and were chosen to bracket 
the effects of friction. The most significant effect of friction, be- 
sides influencing the level of shear stress, is its influence on peak 
stresses. Increasing friction tends to decrease the maximum radial 
stress at G a 0. This effect occurred for all situations studied. 
Friction also increased the maximum circumferential stress near the end 
of the contact region. Friction had some effect on the actual contact 
zone. Because of the effect of friction at G * 0, the cosinusoidal rep- 
resentation (not shown) is closer to the frictionless case than the 
other cases. In fact, with no friction and no clearance, the cosinusoid- 
al assumption is close for this particular plate. 



33 


Table 4 


Location of Maximum Hoop Stress, Contact and 
No-Slip Regions for Various Cases 


Plate 

Fig. 

fixed 

variable 

location of 

contact 

no-slip 


No. 

conditions 

condition 

max a 

8 

arc, B 

arc, a 




pin 






M * 0.2 

ri gi d 

88° 

74° 

< 5° 

A 

6 

6/R * 0.035 

steel 

88° 

74° 

< 5° 



X/R = 0.01 

aluminum 

88° 

74° 

< 5° 




pin 






M = 0.2 

rigid 

74° 

74° 

< 5° 

B 

7 

6/R = 0.035 

steel 

74° 

74° 

< 5° 



X/R = 0.01 

aluminum 

74° 

74 a 

< 5° 




pin 






M = 0.2 

rigid 

74° 

74° 

< 5° 

C 

3 

6/R = 0,035 

steel 

74° 

74° 

< 5° 



X/R * 0.01 

aluminum 

76° 

76° 

< 5° 




pin 






u • 0.2 

rigid 

74° 

74° 

< 5° 

D 

9 

6/R = 0.35 

steel 

74° 

74° 

< 5° 



X/R * 0.01 

aluminum 

76° 

75° 

< 5° 



steel pin 

u 

0 

76° 

71° 


13 

10 

6/R * 0.035 

0.2 

74° 

74° 

< 5° 



X/R = 0.01 

0.4 

77° 

77° 

10 ° 




X/R 






v = 0.2 

0 

86° 

86° 

< 5° 

D 

11 

steel pin 

0.01 

74° 

74° 

< 5° 



6/R * 0.35 

0.02 

69° 

56° 

< 5° 




6/R 






M = 0.2 

0.02 

66° 

61° 

< C° 

0 

12 

steel pin 

0.035 

74° 

74° 

< 5° 



X/R - 0.01 

0.05 

78° 

78° 

< 5° 
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Figure 11 illustrates the effect of pin/hole clearance on the hole- 
edge stresses in plate D. The three levels of clearance illustrated 
are x/R = 0.0, 0.01, 0.02. The value x/R * 0 is often referred to as 
the snug-fit, push-fit, or perfect-fit pin. The value of p used was 0.2 
and 6/R was 0.035. The pin was steel. The most effect pin/hole clear- 
ance has is in the location of the peak circumferential stress. In- 
creasing the clearance moves the peak stress location toward 0 = 0. 

This is a direct consequence of the rapid decrease in contact zone with 
increasing clearance. Increasing the clearance also causes more of a 
region to experience high circumferential stress. Table 4 quantifies 
the effect of clearance on contact angle. As expected, the value of the 
peak radial stress is sensitive to the level of clearance. With less 
contact area the pin naturally loads the contact region more. 

Finally fig. 12 shows the effect of increasing pin displacement on 
the stresses. The pin is steel, the coefficient of friction is 0.2, and 
the clearance is 0.01. Three values of pin displacement are considered: 
6/R * 0.02, 0.035 and 0.05. The most significant effects are with the 
contact arc and the location of the peak circumferential stress. With 
increasing pin displacement the contact arc increases and the location 
of the peak circumferential stress moves toward e = 90°. Obviously in- 
creasing pin displacement Increases the magnitude of the actual stres- 
ses. However, in the nondimensional sense shown in fig. 12, the stress 
magnitudes are not strongly Influenced. The peak nondimensional radial 
stress actually decreases some with Increasing pin displacement due to 
the longer contact arc. 
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Additional Comments 

Two other issues should be addressed before closing. The first is- 
sue deals with the functional behavior of the stresses between colloca- 
tion points. The second issue deals with friction. 

Figures 6-12 were drawn by hand-fairing lines through the numerical 
values of the stresses computed at the collocation points. When a 
Fourier series representation of a function is used in conjunction with 
collocation, there is the issue of behavior of the series representation 
between the collocation points. In addition, here there is the issue of 
the behavior near the end of the contact zone, where the stress state is 
rapidly changing with e. When the numerical values of A k were substitu- 
ted into eq. 7 and N, for example was computed as a function of e, the 
function did indeed oscillate between collocation points. For a typical 
situation, for 6 < @ the amplitude of the oscillation was less than 2% 
of its maximum value. Just beyond 0 3 e the amplitude of the oscilla- 
tion increased, but not markedly. Oscillation to a certain degree was 
expected. Wilson and Goree [24] discussed * dual -series approach to a 
contact problem which resulted in very little oscillation anywhere. 

Such an approach could have been used in this problem. However, nothing 
other than minimal oscillations were experienced with the numerical ap- 
proach here and so the issue of other numerical schemes did not arise. 

Concerning friction, Dundurs and Comninou [25] have shown that at 
the transition point between no-slip and slip, the slope of x r0 , with 
respect to 0, becomes infinite. With the approach here such behavior 
would never be realized and thus the results are in error in that re- 
gard. An alternative approach would have been to build in a special 
friction-induced shear function which did indeed yield an infinite slope 
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of the transition point. However, it is felt that the results presented 
are not in considerable error for not accurately representing this 
effect. 

Concluding Remarks 

This paper has used elasticity solutions and a numerical procedure 
to study the stress distributions around a hole in a pin-loaded ortho- 
tropic plate. In particular, the effects of pin elasticity, friction 
and clearance have been studied. In addition, the effects of the 
plate's elastic properties on the stress distribution has been assessed. 
It can be concluded that, within the context of this study, pin elasti- 
city is not an important variable. Pin/hole clearance is an important 
variable and friction does effect the stress distributions. Further 
studies are warranted in the area of finite-geometry plates, a circum- 
ferentially variable coefficient of friction, or alternatively, a more 
general (non-Coulomb) friction model. Perhaps experimental investiga- 
tions into the slip and no-slip zones are necessary before more general 
friction laws are incorporated. 
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Fig. 1: a) Geometry and coordinate system used in analysis. 

b) Definition of clearance, \, and pin displacement, 6. 
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Fig. 2: Important regions at pin/hole boundary. 
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Fig. 3: Idealization of pin load: a) pin in double shear; b) section 

of pin within plate thickness. 
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Fig. 5: Choosing the correct value of «. 
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Fig. 7: Stresses around hole in plate B. 
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Fig. 8: Stresses around hole In plate C. 
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Fig. 9: Stresses around hole in plate D 
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Fig. 11: Effect of clearance on stresses. 





